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limit is proved to be monotonic (non-monotonic) for FBC y (PBC^). It is 

noted that ground state entropy determinations on infinite strips can be used 

to obtain the central charge for cases with critical ground states. 
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I. INTRODUCTION 



Nonzero ground state entropy, So 7^ 0, is an important subject in statistical mechanics. 
One physical example is provided by ice, for which the residual molar entropy is Sq = 
0.82 ± 0.05 cal/(K-mole), i.e., S /R = 0.41 ± 0.03, where R = N Avog k B [|,|. This is 
equivalent to a ground state degeneracy per site W > 1, since So = ks^W. Such nonzero 
ground state entropy violates the third law of thermodynamics (see, e.g., P,f|). The g-state 
Potts antiferromagnet (AF) exhibits nonzero ground state entropy (without frustration) 
for sufficiently large q on a given lattice A and serves as a useful model for the study of this 
phenomenon. There is an interesting connection with graph theory here, since the zero- 
temperature partition function of the above-mentioned g-state Potts antiferromagnet on a 
graph G satisfies Z(G,q,T = 0)paf = P(G,q), where P(G,q) is the chromatic polynomial 
expressing the number of ways of coloring the vertices of the graph G with q colors such that 
no two adjacent vertices have the same color @. Thus, 

W{[ lim G }, q) = lim P(G, q) 1/n (1.1) 

n— >oo n— >oo 

where n = v(G) is the number of vertices of G. Nontrivial exact solutions for this function 
are known in only a very few cases, all for 2D lattices: the square lattice for q = 3 ||, 
triangular lattice and, for q = 3, the kagome lattice M. Of course, one can use large-g 
series expansions [|T^-|T^], rigorous upper and lower bounds P^3|JT2] -4T4], and Monte Carlo 



simulations (see, e.g., [|16|-|18|| ) . It is also worthwhile to generalize q from Z + to C and study 
q) in the complex q plane for infinite-n limits of various families of graphs, {G} fllDfl - 
|[26|| . On the positive real axis, W({(j}, q) is an analytic function down to a point which we 
denote q c ({G}) |2j. 

Since it is possible to obtain exact analytic solutions for infinitely long strips of 2D 
lattices ||25|| , one has an alternate way to investigate W(A,q) for 2D lattices, namely to 
calculate exactly W on such infinitely long strips of progressively greater widths. These 
W functions for infinitely long strips have interesting analytic structure in their own right, 
which was investigated in detail in Ref. |[25|| . Here we shall use them for a different purpose: 



to investigate how rapidly the 2D thermodynamic limit is approached as the width of the 
strips increases. We find that this approach is quite rapid. Of course, to demonstrate this 
does not require the use of exact analytic results; it can be seen equivalently from numerical 
Monte Carlo measurements on rectangular L x x L y patches after an extrapolation to L x = 00. 
Indeed, Monte Carlo measurements would be the standard method for this purpose since 
they are not limited, as the exact analytic calculations are, to a rather small range of L y 
values. However, the value of discussing this with exact results is that the reader can verify 
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the conclusions directly rather than having to reproduce them with another Monte Carlo 
study. 

The organization of the paper is as follows. In section 2 we discuss some generalities of our 
approach. In section 3 we prove that if one uses free (periodic) boundary conditions in the 
y direction transverse to the length of the infinite strip, then W approaches its 2D thermo- 
dynamic limit monotonically (nonmonotonically) . Section 4 contains the numerical results 
for strips of the square, triangular, and honeycomb lattices with free transverse boundary 
conditions, while section 5 contains analogous results for periodic transverse boundary con- 
ditions. In section 6 we remark on how these strip studies can be used to determine the 
central charge for cases with critical ground states. Our conclusions are presented in section 
7. 



II. W ON STRIP GRAPHS AND THE APPROACH TO 2D THERMODYNAMIC 

LIMIT 

The usual thermodynamic limit of the Potts antiferromagnet or other statistical mechan- 
ical model on the 2D lattice A involves taking L x — > oo and L y —>■ oo with fixed L y /L x = p, 
where (p 7^ 0,oo). The question of how various thermodynamic quantities approach their 
2D limits as a function of p has been of interest for many years (e.g., Ref. [27] for the Ising 



model). As noted in the introduction, a different way to approach the 2D thermodynamic 
limit is via a sequence of infinitely long strips of progressively greater and greater widths. 
That this is different is clear from the fact that for each such strip, regardless of how large 
L y is, the ratio L y /L x = 0. We picture the strip graphs as extending longitudinally in the 
horizontal (x) direction and having a width of L y vertices in the vertical direction. A priori, 
it is not clear that this different approach will yield results that are useful to the study of 
the 2D thermodynamic limit, because, for a given thermodynamic quantity of interest, these 
results might be dominated by the fundamentally ID nature of the infinite strip. Indeed, 
to illustrate a case where it is not useful, consider a model, such as a discrete ferromagnet, 
which has a second-order phase transition at some critical temperature T C (A) on a 2D lattice 
A, and assume that there is no exact solution of this model. If one were to try to employ 
exact solutions of the model on infinitely long strips of lattice type A to determine T C (A) 
for the 2D lattice, one would get the ID result T c = for any finite value of L y . Hence, for 
T C (A) this method would not give any useful information. However, as we shall show, the 
situation is very different with the ground state entropy So (A); for this quantity, one can use 
results on infinite strips to get quite accurate values even for rather modest strip widths. 
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The existence of the thermodynamic limit for the 2D lattice A means that the maximal 
finite real q where W(A,q) is nonanalytic, <7 C (A), is independent of the boundary conditions 
used in taking the 2D thermodynamic limit ||22|| . Let us denote the W function for the 
L x x L y strip of the lattice of type A as W(A(L X x L y ), BC X , BC y ,q). We observe here 
that for physical (positive integral) q > q c (A), in the limit L x — > oo, this W function is 
independent of the boundary conditions used in the x direction. This is also true for real 
q > g c (A) (and more generally, in the region of the complex q plane denoted Ri in our 
previous studies p2|). We thus introduce a more compact notation for the W function on 
infinitely long strips: 

W(A(L y ), BCy, q) = lim W(A(L X x L y ), BC X , BC y , q) (2.1) 

Indeed, let A^ be an infinite d-dimensional lattice and A^_ 1 ^ d be a slab of a (d — 1)- 
dimensional lattice, infinite in d — 1 dimensions and of finite thickness L d in the cPth di- 
mension. For physical q > q c (A d ), the value of W(A d _ 1Ld , q) is independent of the boundary 
conditions used for the (d — 1) directions when taking Lj — > oo for 1 < j < d — 1. As we 
have discussed before j22j| , this is not true for all q G C; however, here we deal only with 
physical q values. 

III. ISSUE OF MONOTONICITY OF APPROACH TO 2D THERMODYNAMIC 

LIMIT 

In this section we show that for free (periodic) boundary conditions in the y direction, W 
for infinitely long strips of width L y approaches its 2D thermodynamic limit monotonically 
(nonmonotonically) as L y — ► oo. 

A. Monotonic Approach for FBC y 

We begin with the case of free boundary conditions and state the following theorem: 
Theorem 1 

Let Ad-i,L d denote a regular lattice graph of infinite extent in d — 1 dimensions and width 
(thickness) L d in the cfth dimension. Let the boundary conditions in the cfth direction be 
free and the boundary conditions in each of the first d — 1 be (separately) free or periodic. 
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(Note that A d -i,i = A d _i and A rf _i j00 = A d .) Then for fixed q > q c (A d ), W (A d _ hLd , q) is a 
monotonically decreasing function of L d for 1 < L d < oo. 

Proof 

We shall prove the theorem for the case d = 2; its generalization to d > 3 will be obvious. 
We consider a finite strip graph of the lattice, of size L x x L y , where the longitudinal direction 
is x. Assume that one has made an allowed coloring of this graph. Now connect another layer 
of sites to the layer that formerly constituted the top layer of sites on the strip. The coloring 
of this new layer of sites imposes additional constraints on the coloring of the original strip, 
and excludes a certain subset of what were previously allowed colorings. Thus, the fraction of 
sites on the augmented graph that have more constraints increases; i.e., the sites on the upper 
and lower edges, which have fewer constraints on their coloring because of the free transverse 
boundary conditions, constitute a progressively smaller fraction of the total number of sites 
as L y increases. Hence, the chromatic polynomial per site, P(A, [L x x L y ],g) 1//n , decreases. 
This inequality holds for arbitrary L x . Taking the limit as L x — > oo and using the definition 
( |1.1| ), one obtains the theorem for the case d = 2. A straightforward generalization of this 
argument proves the theorem for d > 3. □ 

A corollary of this theorem is that if one compares W on two infinite lattices of the same 
type and of different dimensions, such as d- dimensional cartesian lattices then, for fixed 

q > <? C (A), 

W(A d ,q) < W(A d/ ,q) if d>d! (3.1) 

To prove this corollary, one starts with d' = d — 1 and (i) constructs A^ from A^-i by impos- 
ing free boundary conditions in the d'th direction and adding layers of vertices in this cfth 
direction, (ii) uses the monotonicity relation of theorem 1 for the quantities W(A d _i^ d , q), 
and (iii) takes the number of added layers in the (i'th direction to infinity to get A^. The 
monotonicity relation for the infinite lattices (|3.1|) was previously noted by Chow and Wu 
Q. It is important to observe that the monotonicity relation (|3.1| ) does not imply our mono- 
tonicity theorem 1. This is clear from the fact that the inequality ( [3.1|) holds independent 
of the boundary conditions that one uses to define the respective thermodynamic limits on 
A d and A^-i, whereas, on the contrary, the inequality in our theorem 1 does not apply if one 
uses periodic boundary conditions for the (i'th direction of the (d — l)-dimensional strip or 
slab of width L d (see below). 
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B. Non-monotonic Approach for PBC,, 

Next, we show that a similar monotonicity result does not hold if one imposes periodic 
boundary conditions in the cfth direction. This is clear from the proof, since the greater 
freedom in coloring the sites on the boundary layer in the cfth dimension played a crucial 
role, but if one imposes periodic boundary conditions in the cfth direction, there is no such 



boundary layer. The simplest illustration is provided by the case d = 1, for which p2 
g c (Ai) = 2 and 

W(A 1 ,q) = q-l (3.2) 

For free boundary conditions, the function that enters on the right-hand side of eq. ( |1 . 1|) is 

P((A!) n , FBC, q) 1 ^ = q l / n {q - l) 1 "* (3.3) 

For fixed q > q c (Ai), this is a monotonically decreasing function of n as n increases from 1 
to infinity. However, if we impose periodic boundary conditions, i.e. deal with an n-vertex 
circuit graph C n [], then, the function that enters on the right-hand side of eq. (|1.1|) is 

r i V n 

P((A 1 ) n , PBC, q) l ' n = (q - 1) [l + {-l) n {q - l)"^ 1 )] (3.4) 

This is a non-monotonic function of n. For example, for the lowest value of q where the 
ID Potts AF has nonzero ground state entropy, viz., q = 3, for which the n — > oo limit is 
W(Ai,q = 3) = 2, eq. fl3.4|) exhibits the non-monotonic behavior indicated by the values 
6 1 / 3 = 1.817.. for n = 3, (18) 1 / 4 = 2.060.. for n = 4, (30) 1 / 5 = 1.974.. for n = 5, (66) 1 / 6 = 
2.010.. for n = 6, etc. Similar non-monotonic behavior occurs for higher values of q. Looking 
at subsequences, we find that P((Ai) n , PBC, g) 1//n is a monotonically increasing function 
of n for odd n > 3 and a monotonically decreasing function of n for even n > 2. This 
is connected with the fact that the circuit graph [(Ax) n , PBC] = C n with odd (even) n 
has chromatic number % = 3 (x = 2). The different behaviors of P((Ai) n , FBC, q) 1 ^ and 
P((Ai) n , PBC, g) 1//n can be seen in a more general context by analytically continuing eqs. 



1 Parenthetically, we note that C n is only a (proper) graph for n > 3 since the strict mathematical 
definition of a graph forbids (i) any multiple bond connecting a given pair of vertices (present for 
C n= 2) and (ii) any bond going out from a given vertex and looping back to the same vertex (present 
for C n= \). This is not important for our demonstration of non-monotonicity. 
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( p. 3D and (|3.4Q from n G Z + to n G R+ and plotting them as functions of n (in the second 
case, since P(C n ,q) — (q — 1)" + (— l) n (q — 1) is complex for n Z, we plot \P(C n , q)\ 1 ' n ). 
This is shown in Fig. |l[ One notices that although eq. ( |3.4|) for periodic boundary conditions 
behaves non-monotonically, it approaches the n = oo value W(Ai,q = 3) = 2 considerably 



more rapidly than the FBC expression, eq. (|3.3| ). As one increases q beyond 3, the first 
maximum in \P(C n , g)| 1//ra moves slightly leftward, and the oscillations damp out faster. 
As we shall show in the tables below, a similar difference holds between the behavior of 
W(A(L y ), BC y , q) for FBC^ and PBC r 




10.0 



FIG. 1. Plots of (a) eq. (3^3) and (b) eq. (3^4) for q — 3, as functions of real n, together with physical 
values for integer n. Horizontal line is the asymptote, W(Ai, q = 3) = 2. 



IV. QUANTITATIVE RESULTS FOR STRIPS WITH FBC y 

One would like to go beyond the general inequality in Theorem 1 to obtain an explicit 
numerical determination of the dependence of W on L^. We do this here for d = 2 and, 
in particular, for the square (sq), triangular (t), and honeycomb (he) lattices. For the strip 
graph of each type G s , we define the ratio 

Rw(HLl) , BCl ,^^^ (4.1) 

In the appendix we list the exact analytic expressions for W(A(L y ), FBC y , q) for A = sq, t, he 
and the L y values used here. In Table | we show a numerical comparison for strips of the 
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square lattice (along the row direction) for 1 < L y < 4 and 3 < q < 10. The exact value 
W(sq,q = 3) = (4/3) 3 / 2 = 1.53960... is from Ref. @, while the values of W{sq, q) for 
4 < q < 10 are from our Monte Carlo measurements in Ref. [ 22| . Using the conservatively 
quoted uncertainties that we gave for the Monte Carlo measurements, it would follow that 
the corresponding uncertainties in the ratios (|4.1| ) are ~ (3 — 4) x 10 -4 ; with less conservative 
estimates of uncertainties in the Monte Carlo measurements, the resultant uncertainties in 
these ratios would be smaller. In Tables || and |T| we give the analogous comparisons for 
strips of the triangular lattice of widths L y = 2, 3, 4 and of the honeycomb lattice for widths 
L y = 2 and 3. In all of these cases, one observes that, for fixed q, the agreement with the 
infinite-lattice value gets better as the width increases and for fixed width, the agreement 
gets better as the value of q increases. These comparisons show that the approach to the 2D 
thermodynamic limit is reasonably rapid even for free boundary conditions in the transverse 
direction. For example, for an oo x 4 strip of the square lattice for 3 < q < 5, the W values 
are within about 5 % of their respective values for the infinite 2D lattice. 



7 



TABLE I. Comparison of values of W(sq(L y ), FBC yi q) with W(sq, q) for 3 < q < 10. For each value of 
q, the quantities in the upper line are identified at the top and the quantities in the lower line are the values 
of Rw(sq(L y ) 7 FBCy,q). The FBC y is symbolized as F in the table. 



q 


^(1)^,9) 


W(sq(2),F,q) 


W(sq(3),F,q) 


W(sq(4),F,q] 


) W(sq,q) 


W(sq,q) s 


W(sq,ql 


3 


2 


1.73205 


1.65846 


1.624945 


1.53960.. 


1.53960.. 


1.50000 




1.299 


1.125 


1.077 


1.055 


1 


1 




4 


3 


2.64575 


2.53800 


2.48590 


2.3370(7) 


2.3361 


2.33333 




1.284 


1.132 


1.086 


1.064 


1 


~ 


~ 


5 


4 


3.60555 


3.48304 


3.42336 


3.2510(10) 


3.2504 


3.25000 




1.230 


1.109 


1.071 


1.053 


1 


_ 


~ 


6 


5 


4.58258 


4.45136 


4.38717 


4.2003(12) 


4.2001 


4.20000 




1.190 


1.091 


1.060 


1.0445 


1 


~ 


~ 


7 


6 


5.56776 


5.43073 


5.36348 


5.1669(15) 


5.1667 


5.16667 




l.lOl 


1.0 (0 


1.051 


l.Ooo 


1 






8 


7 


6.55744 


6.41623 


6.34677 


6.1431(20) 


6.1429 


6.14286 




1.1395 


1.067 


1.0445 


1.033 


1 






9 


8 


7.54983 


7.40548 


7.33434 


7.1254(22) 


7.1250 


7.12500 




1.123 


1.060 


1.039 


1.029 


1 






10 


9 


8.54400 


8.39720 


8.324745 


8.1122(25) 


8.1111 


8.11111 




1.109 


1.053 


1.035 


1.026 


1 
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TABLE II. Comparison of values of W(t(L y ), FBC y ,q) with W(t, q) for 3 < q < 10. For each value of 
q, the quantities in the upper line are identified at the top and the quantities in the lower line are the values 
of Rw{t{L y ), FBC y ,q). The FBC y is symbolized as F in the table. 



q 


W(t(2),F,q) 


W(t(S),F,q) 


W(t(4),F,q) 


W(t,q) 


W(t,q) t 


4 


2 


1.77173 


1.67619 


1.46100 


1.333333 




1.369 


1.213 


1.147 


1 


_ 


5 


3 


2.72998 


2.60495 


2.26411 


2.250000 




1.325 


1.206 


1.151 


1 


_ 


6 


4 


3.71457 


3.579715 


3.20388 


3.200000 




1.248 


1.159 


1.117 


1 


_ 


7 


r. 



A 70^71 
l ±. 1 UO ( 1 


4t. OOOIO 


A 1 RSI Q 


^t. IDODD i 




1.200 


1.129 


1.095 


1 




8 


6 


5.69974 


5.55530 


5.14358 


5.142857 




1.167 


1.108 


1.080 


1 




9 


7 


6.695395 


6.54810 


6.12539 


6.125000 




1.143 


1.093 


1.069 


1 




10 


8 


7.69208 


7.54259 


7.11134 


7.111111 




1.125 


1.082 


1.061 


1 





TABLE III. Comparison of values of W(hc(L y ),FBC y ,q) with W(hc, q) for 3 < q < 10. For each value 
of q, the quantities in the upper line are identified at the top and the quantities in the lower line are the 
values of Rw(hc(L y ),FBC y ,q). The FBC y is symbolized as F in the table. 



q 


W(hc(2),F,q) 


W(hc(S),F,q) 


W(hc,q) 


W(hc, q) ser . 


W(hc,q) e 


3 


1.82116 


1.76567 


1.6600(5) 


1.6600 


1.658312 




1.097 


1.064 


1 


_ 


_ 


4 


2.79468 


2.72942 


2.6038(7) 


2.6034 


2.603417 




1.073 


1.048 


1 


_ 


_ 


5 


3.78389 


3.71448 


3.5796(10) 


3.5795 


3.579455 




1.057 


1.038 


1 


_ 


_ 


6 


4.77760 


4.70568 


4.5654(15) 


4.5651 


4.565085 




1.046 


1.031 


1 


_ 


_ 


7 


5.77336 


5.69973 


5.5556(17) 


5.5553 


5.555278 




1.039 


1.026 


1 






8 


6.77028 


6.69539 


6.5479(20) 


6.5481 


6.548095 




1.034 


1.023 


1 






9 


7.76793 


7.69208 


7.5424(22) 


7.5426 


7.542587 




1.030 


1.020 


1 






10 


8.76607 


8.68945 


8.5386(25) 


8.5382 


8.538222 




1.027 


1.018 


1 







V. QUANTITATIVE RESULTS FOR STRIPS WITH PBCy 

In Tables [IV] and |V| we present similar results for infinite strips with periodic boundary 
conditions in the transverse (y) direction. The exact analytic expressions that we use for 
these tables are given in the Appendix. As is mentioned in the Appendix, for a strip of the 
square lattice with PBC^ and cross sections forming triangles, depending on one's labelling 
conventions, this corresponds to L y = 3 or L y = 4, where in the latter case, one interprets 
the periodic boundary conditions as identifying the top and bottom vertices for each value 
of x. A similar comment applies for a strip with PBC y and transverse cross sections forming 
squares. For the table, we use the convention of choosing the smaller of the respective values 
of L y . We find that for a given q, W approaches its 2D value W(A, q) much more rapidly 
with periodic rather than free transverse boundary conditions: for the modest width L y = 4, 
W is within O(10~ 3 ) of its 2D value for moderate q. The finding that the periodic boundary 
conditions in the transverse direction yield a more rapid approach to the 2D thermodynamic 
limit than the free boundary conditions is not, in itself, a surprise; this is in accord with a 
wealth of past experience with statistical mechanical models on finite-size lattices. What is 
remarkable is how rapid in absolute terms this approach is. Of course, one can also consider 
larger values of L y , but the strikingly rapid approach to the 2D thermodynamic limit is 
already fully demonstrated by the range of L y that we have considered. 



TABLE IV. Comparison of values of W(sq(L y ), PBC y , q) with W(sq, q) for 3 < q < 10. For each value 
of q, the quantities in the upper line are identified at the top and the quantities in the lower line are the 
values of Rw(sq{L y ), PBC yi q). The PBCy is symbolized as P in the table. 



Q 


W(sq(3),P,q) 


W(sq(A),P,q) 


W(sq,P,q) 


3 


1.25992 


1.58882 


1.53960.. 




0.8183 


1.032 


1 


4 


2.22398 


2.37276 


2.3370(7) 




0.9516 


1.015 


1 


5 


3.17480 


3.26878 


3.2510(10) 




0.9766 


1.0055 


1 


6 


4.14082 


4.21082 


4.2003(12) 




0.9858 


1.002505 


1 



11 



10 



5.11723 
0.9904 
6.10017 
0.9930 
7.08734 
0.9947 
8.07737 
0.9957 



5.17377 

1.0013 
6.14792 

1.0008 
7.12881 

1.0005 
8.11409 

1.0002 



5.1669(15) 
1 

6.1431(20) 
1 

7.1254(22) 
1 

8.1122(25) 
1 
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TABLE V. Comparison of values of W{t{L y ), PBC y , q) with W(t, q) for 4 < q < 10. For each value of 
q, the quantities in the upper line are identified at the top and the quantities in the lower line are the values 
of Rw(t(L y ),PBCy,q). The PBCj, is symbolized as P in the table. 



q 


W(t(S),P,q) 


W(t(4),P,q) 


W(t,q) 


4 


1.58740 


1.18921 


1.46100 




1.0865 


0.8140 


1 


5 


2.35133 


2.21336 


2.26411 




1.0385 


0.9776 


1 


6 


3.23961 


3.185055 


3.20388 




1.0112 


0.9941 


1 


7 


4.17934 


4.15965 


4.16819 




1.0027 


0.99795 


1 


8 


5.14256 


5.13936 


5.14358 




0.99980 


0.9992 


1 


9 


6.11803 


6.12324 


6.12539 




0.99880 


0.99965 


1 


10 


7.10059 


7.11027 


7.11134 




0.99849 


0.99985 


1 



VI. CASES WITH CRITICAL GROUND STATES 



For certain 2D lattices A and values of g, the g-state Potts antiferromagnet has a critical 
ground state, i.e., as T — > 0, a correlation length £ defined, say, by a spin-spin correlation 
function, goes to infinity. Normally, in statistical mechanics, for a given dimensionality d and 
symmetry group G, second-order phase transitions can be described by a universality class 
representing a fixed point of the renormalization group. Conformal field theory methods 
have provided a powerful way to understand these universality classes and the associated 
critical exponents in terms of Virasoro algebras with given central charges and scaling di- 
mensions |f28fl . In addition to phase transitions involving ferromagnetic long range order at 
low temperatures, this is also true of antiferromagnetic transitions on bipartite lattices, but 
the situation is more complicated on nonbipartite lattices, as is illustrated by the fact that 
the isotropic Ising antiferromagnet on the triangular lattice has no finite-temperature phase 
transition but is critical at T = 0. 

The q = 3 Potts antiferromagnet on the square lattice has a critical ground state with 
central charge c = 1, as a consequence of the fact that at T = this model can be mapped 
to a critical six- vertex model ||. From the exact solution in Ref. |§, it can be argued that 
the q = 4 Potts antiferromagnet on the triangular lattice is also critical, which is closely 
related to the fact that the q = 3 Potts antiferromagnet on the kagome lattice is critical at 



T = ||29|| . We recall that given the Virasoro algebra with central extension 

[L m , L n ] = (m - n)L m+n + —m{m 2 - l)5 m+nfi (6.1) 
and the corresponding Kac-Moody algebra realized at level k 



[■C J b n\ = C abC J C m + n + TT^Wo (6.2) 



with structure constants c a bc, as connected via the Sugawara relation (e.g., p8|j ) 

^ oo 



it follows that 



m=— oo 



dim(g) 



(6.4) 



C 2 (g)/k + l 

where 02(g) is the quadratic Casimir operator for the algebra g. In particular, 

g = su(M) k=1 c = M - 1 (6.5) 
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Hence, from eq. (|6.5| ) together with the finding [29[] that the Kac-Moody algebra is su(3)k 



=1 



for the T = q = 3 Potts AF on the kagome lattice [2J|, it follows that c = 2 for this critical 
ground state. Given that there is a close connection between the Potts antiferromagnets 
with q = 3 on the kagome lattice and with q = 4 on the triangular lattice, which leads to 
the relation W(kag, q — 3) = W(tri, q = 4) 1 / 3 ||, this suggests that this value of c = 2 also 
holds for the T = 0, q = 4 Potts AF on the triangular lattice. 

Here we would like to point out that determinations of the ground state entropy on 
infinitely long strips of finite width can be used to obtain the central charge c for Potts 
antiferromagnets with critical ground states. If one considers the model on a lattice of size 



L x x L y , then, in the limit as L x — > oo, one has [31 



fstri P ,L y — fbulk + + yrr + 0(L y 3 ) (6.6) 

where f sur f. = is nonzero (zero) for free (periodic) boundary conditions in the y direction 
and 

I c for PBCL 

A = { 6 V (6.7) 

for YBC y 

For the critical ground states of interest here, viz., q = 3, 4, 3 on the square, triangular, and 
kagome lattices, respectively, as well as other possible 2D cases, the Potts antiferromagnet 
exhibits ground state entropy without frustration, and the reduced free energy (per site) / = 
liniAr^oo iV _1 In Z is given simply by the ground state entropy: /(A, q) paf = So(A, q) pAF/ks, 
Hence, eq. (|6.6j) becomes 



S s trip,L y — Sbulk H + Y2 ^(-^2/ 3 ) 



Thus, calculations of S st rip,L y for several different values of L y can yield c. Normally, one 
would do this via the most general and robust method, namely Monte Carlo simulations. For 
small L y values, we have shown that it is actually possible to get exact analytic results, but 
this method is not competitive with Monte Carlo simulations for strips with larger values 
of L y . One might note in passing that for the q = 3 Potts AF on the infinite square strip 
with PBCy and cross sections forming squares, taking L y = 4, using the fact that S sur f, = 

-3 



in this case, and dropping terms of order L y in eq. (|6.8| ), we obtain the estimate c = 0.96, 
quite close to the exact value c = 1. 
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VII. CONCLUSIONS 



In summary, we have studied a different type of approach to the 2D thermodynamic 
limit for the ground state entropy, or equivalently, the ground state degeneracy per site, 
W, of the g-state Potts antiferromagnet, using infinitely long strips of increasing widths. 
We have found that the approach of W to its 2D thermodynamic limit is quite rapid; for 
moderate values of q and widths L y ~ 4, W(AL y ,q) is within about 5 % and 0(1CT 3 ) of 
the 2D value for free and periodic boundary conditions, respectively. We have also proved 
that the approach of W to the 2D thermodynamic limit is monotonic (non-monotonic) for 
free (periodic) boundary conditions in the transverse direction. Finally, we have noted that 
these ground state entropy determinations on infinite strips can be used to obtain the central 
charge for cases with critical ground states. 



We are grateful to Prof. M. Rocek for the collaborative work on Ref. [25] . This research 
was supported in part by the NSF grant PHY-97-22101. 



VIII. APPENDIX 

We gather together here the exact analytic formulas on which our numerical tables are 
based. It should be emphasized that the entries in these tables and the resultant conclusions 
about the rapidity of the approach of W to the 2D thermodynamic limit for infinitely long 
strips with free or periodic transverse boundary conditions could also have been obtained 
using purely numerical Monte Carlo calculations. The usefulness of the analytic formulas 
(which are elementary for L y = 1,2) is just that they enable one to check the results more 
directly. 



A. Square Lattice, FBC^ 

For infinitely long strips of the square lattice with FBC y , we have 

W(sq(L v = l),FBC v ,q) = q-l (8.1) 

W(sq(L y = 2), FBC y , q) = (q 2 - 3q + 3) 1 / 2 (8.2) 



16 



W(sq(L y = 3),FBC y ,q) = 2- 1 / 3 



(g-2)(g 2 -3g + 5) 



+ 



(g 2 - 5q + 7)(g 4 - 5g 3 + llq 2 - Ylq + 



n 1/2 



1/3 



W^sq^Ly = 4), FBCy, q) is given by the maximal root of the cubic equation 

£ 3 + frs 9 (4),l£ 2 + &s 9 (4),2^ + &sg(4),3 = 

where the coefficients 6 a g(4),jfcj fc = 1, 2, 3 were listed in Ref. pjfl . 



B. Triangular Lattice, FBC 4 

For the triangular lattice strips with FBC^, we have 

W{t{L y = 2),FBC y ,q)=q-2 



W(t(L y = 3),FBC y ,q) = 2- 1 / 3 



(q 3 - 7q 2 + 18q - 17) 



+ 



q 6 - Uq 5 + 81q 4 - 250g 3 + 442g 2 - 436q + 193 



1/2' 



n l/3 



W(t(L y = 4), FBCy], q) is given by the maximal root of the quartic equation 

^ + ^(4),l£ 3 + ^(4), 2 e 2 + &t(4), 3 f + &t(4),4 = 

where the btuvk, k — 1, ..,4 were listed in Ref. fl2"5f . 



C. Honeycomb Lattice, FBC y 

For the honeycomb lattice strips with FBC^, we have 

W(hc(L y = 2), FBCy, q) = (q 4 - 5q 3 + Wq 2 - lOq + 5) 1/4 
WihciLy = 3), FBCy, q) is given by the maximal root of the cubic equation 

£ 3 + fr/ic(3),l£ 2 + hc(3),2^ + bhc(3),3 = 

where the b hc( 3) tk , 



k = 1, 2, 3 were listed in [25 
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D. Square Lattice, PBC^ 

We first consider a strip of the square lattice with PBC^ and transverse cross sections 
forming triangles. Depending on one's labelling conventions, this corresponds to L y = 3 
or L y = 4, where in the latter case, one interprets the periodic boundary conditions as 
identifying the top and bottom vertices for each value of x. We calculate 

W{sq{L y = 3), PBC y , q) = (q 3 - 6q 2 + Uq - 13) 1/3 (8.10) 



For the next larger size, i.e. transverse cross sections forming squares, corresponding to 
Ly = 4 or Ly = 5 in the respective labelling conventions described above, the W function is 
given by f25f 



W(sq(L y = 4), PBCy, q = 2" 1 / 4 



(q 4 - 8q 3 + 29q 2 - 55q + 46) 



+ 



q 8 - 16q 7 + 118<f - 526g 5 + 1569g 4 - 3250g 3 + 4617g 2 - 4136g + 1776 



1/2' 



■,1/4 



(8.11) 



E. Triangular Lattice, PBC, 



We next consider a strip of the triangular lattice with PBC^, represented as a square 
lattice with additional diagonal bonds from, say, the upper left to lower right vertices of 
each square. For the case where the transverse cross sections form triangles, corresponding 
to L y = 3 or L y = 4 in the above labelling conventions, we calculate 



W(t(L y = 3), PBCy, q) = (q 3 - 9q 2 + 29q - 32) 1/3 
For the next larger size, with transverse cross sections forming squares, W is PS - 



(8.12) 



W(t(L y = 4), FBCy], q) = 2~ 1 / 4 (g - 3) 1 



(q 3 - 9q 2 + 33q - 4£ 



+(?-4) 



q 4 - 10q 3 + 43q 2 - 106g + 129 



1/2' 



-, 1/4 



(8.13) 
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